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Abstract. A first order phase transition is found in a model which was introduced originally 
by Murthy and Shankar [Phys. Rev. B 60, 6517 (1999)] to describe systems of generalised exclusion 
statistics. I characterise the phase transition in the canonical and grandcanonical ensebles for the 
case when the statistical exclusion parameter is 1, which corresponds to the Fermi exclusion statistics. 
We observe that in the grandcanonical ensemble the phase transition has no latent heat, but it has a 
finite jump in the particle number. In canonical conditions-when the particle number is held fix— the 
internal energy is discontinuous at the transition. 
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1 INTRODUCTION 

The exclusion statistics that it is intermediate between the Bose and the Fermi statis- 
tics, called fractional exclusion statistics (FES), was introduced by Haldane in Ref. [T] 
to describe, among other systems, quasiparticle excitations in fractional quantum Hall 
effect and spinon excitations in spin-^ quantum ferromagnets. For this he assumed 
that the many-particle system has a finite-dimensional Hilbert space whos dimension 
varies linearly with the particle number in the system. 

The thermodynamic properties of systems with FES in the thermodynamic limit 
were clarified mainly by Isakov [2] and Wu [3]. 

One of the interacting models that leads to FES was introduced first by Murthy 
and Shankar in Ref. [4] and then used with slight variations in different contexts by 
various authors (see for example [6l [3 [8] ) . To explain the model, let us assume 
for the beginning that we have an ideal gas of bosons, with the single particle levels 
denoted by e,-, i — 0, 1, . . . and e; < ej+i for any i. I also assume, for the sake of 
clarity, that the energy levels are equidistant and I introduce the density of states, 
a = (fii+i — e-i) -1 - The interaction between particles changes the energy levels in such 
a way that the energy level i changes its energy to 

ii = e(e») = ej +a~ 1 gN i - + a^ho max(rij —1,0), (1) 

where rij is the numer of particles on the level i, while 7Vi_ = X^<i n j ^ s the number 
of particles below the level i. By max(a, b) I denote the maximum of a and b. The 
constants g and Iiq are adimensional interaction parameters. In order to ensure the 
stability of the system, I impose from the beginning ho > and g > 0. Obviously, 
ho — g = corresponds to the ideal Bose gas. 

Let us now analyse the (quasi)continous limit. In this limit we can relax the 
condition of having an equidistant single particle spectrum and describe the system 



only by the average density of states a. Then for any e > 0, if there are no macro- 
scopically populated energy levels-i.e. in the absence of condensation-, n(e) is much 
smaller than iV c _ and we can write 

e(e) = e + g f n(e)de . (2) 



Jo 

Even for e = Owe can neglect the microscopic offset a hoii(e = 0) of the groundstatc- 
which has no contribtion to the thermodynamics of the system-and write e(0) = 0. 

Now we can prove the generalised exclusion statistics character of the system. 
For this, I will take the arbitrary interval [£\, e u ] , which contains a large enough number 
of particles, nrg 1 g u ], and calculate the number of states it contains. The energies ei 
and e u are related to the energy levels in the noninteracting system, ei and e u , by the 
relations 

ei = e\ + a~ 1 \gNi l _ + h o max(n Sv 0)] , (3) 
e u = e u + a~ 1 [g(N il _ + n[g b g u ]) + /i max(n ?a , 0)] . (4) 

For enough many energy states in the interval (i^gj 3> ng v n^) I can write the 
number of states between ei and e u as 

d[e lt £ u ] = cr(e u - q) = a(e u - q) - n^^g . (5) 

Relation ^ is linear in the number of particles in the interval and gives an effective 
exclusion statistics parameter a = g 5, 6, 7, 8J. The case a = corresponds to Bosons, 
whereas a = 1 corresponds to Fermions. Therefore the interacting gas is described as 
an "ideal" gas of particles of generalised statistics, with the same, constant density 
of states a. 

The thermodynamics of such a system is obtained by splitting the energy axis 
into intervals which are small enough, so that the energy variations caused by different 
rearrangements of particles in any of the intervals are insignificant, but which are also 
sufficiently large, so that each interval contains a large number of particles and energy 
levels. If we number these energy intervals by capital letters (e.g. / = 0, 1, . . .) and the 
interval / contains Sni particles and di energy levels, then the number of microstates 
obtained by rearranging the particles in this interval is 

(5m + di)\ 

Wi = 5nMA (6) 

Summing up the number of microstates, multiplied by their statistical weights, one 
obtains the partition function, from which, by taking the logarithm and using the 
Stirling formula, one calculates the thermodynamic potential and entropy (see for 
example 13 IH]) ■ The thermodynamic potentials have been put into an universal 
form for any a by writing them in terms of polylogarithm functions [S] , as previously 
done independently for Bosons (a — 0) and Fermions (a = 1) by Lee [10] and Viefers, 
Ravndal, and Haugset[ll . For example the grandcanonical potential, tt, and the 
entropy, S, can be written in the form 

n = -U= 1 —^— + (k B T) 2 aLt 2 (-y Q ), (7a) 
z a 

S = -klTa[2Li 2 (-y )+\og(l + y )logy }, (7b) 
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if we define y by the equation (1 + yo) 1 ~ a /yo = e~^l kliT . In Eq. (fTaj) TV is the total 
particle number in the system and /x is the chemical potential. 

Equation (|7b[) shows that the entropy of any system of constant DOS is indepen- 
dent of the exlcusion statistics parameter a [S]. From here follows that in canonical 
conditions the systems of any exclusion statistics, but of the same, constant density of 
states, are identical from the thermodynamic point of view. This property has been 
given the name of thermodynamic equivalence of systems of constant DOS [TU] . 

1.1 THERMODYNAMIC EQUIVALENCE 

Apparently ho plays no role in the model since does not enter any of the macroscopic 
equations ([5|), (JH|), or (J7|). For this reason it has been in general neglected [7], if not 
entirely discarded [H O O O [12] . After all, a thermodynamic Bose system is usually 
thought of as an infinite Bose gas, with a — > oo, N — > oc, and finite density p = N/a. 
Such a gas does not condense at any finite temperature, so there is no reason to 
consider h Q . 

Moreover, the model can be built from a fermionic perspective. In [5] I showed 
that the microscopic reason for the thermodynamic equivalence between a Bose and a 
Fermi system of the same, constant density of states is the possibility to make a one- 
to-one mapping between the microscopic configurations of bosons and the microscopic 
configurations of fcrmions, with the same excitation energy |13j . By this mapping, 
the bosons that lay on the same energy level are transformed into a goup of "close 
packed" fermions, i.e. into a group of fcrmions which occupy all the states into a 
single-particle energy interval (more details are given in Refs. H3|). In such a 
case, the single particle energies of the interacting Fermi gas corresponding to Eq. (TTJ) 
should be written as 

ei = e( Ci ) = e, + tr-y ^JV,, + a-'h^P max^ -1,0), (8) 

with g( F ) = 1 — g and h,Q — ho — 0.5. In this description, rij is the number of 
fermions in the "close-packed" group that contains the particle on the level i. Since 
at temperature T in any small energy interval <5e, centered at e, there are on average 
aSe/{l + exp[/3(^i — e)]} 0) holes, we cannot expect that a macroscopical interval 
along the energy axis will be "close-packed," so from this perspective also-maybe even 
more clearly-the parameter ho does not appear to have any macroscopic significance. 

1.2 THE ROLE OF h 

Infinite systems of constant density of states are indeed not condensed at any finite 
temperature, as stated before. Nevertheless, in finite systems, a significant fraction 
of particles may accumulate on the ground state at low enough temperature in Bose 
systems, or it can arrange in a close-packed (or degenerate) configuration on the lowest 
energy levels in Fermi systems |14) . Therefore, at low enough temperature we should 
take into account separately the population of the ground state in Bose systems or 
the corresponding degenerate fermionic subsystem and I shall refer to them as the 
Bose or the Fermi condensates [14] . 

By doing so, in Ref. [H] I showed that if ho < 0.5 the system undergoes a 
first order phase transition for any a. In this paper I will describe in detail the 
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phase transition for a = 1. In Sectional calculate the partition functions and I 
show how the phase transition occures in grandcanonical (Section 12. 1[) and canonical 
systems (Section 12. 2[) . In Section [3] I calculate the discontinuities that appear in 
the thermodynamical functions at the phase transitions. I shall show that in the 
grandcanonical ensemble the internal energy is continuous at the transition, while the 
particle number has a jump. However, in the canonical ensemble, since the particle 
number is conserved, the internal energy is discontinuous at the phase transition. 
Last section is reserved for conclusions. 

2 THE PHASE TRANSITION 

2.1 THE GRANDCANONICAL ENSEMBLE 

From now on I shall discuss the case a = g = 1. Although the analysis is equally 
easy from Bose and Fermi perspectives, for easier reference to the calculations in [5] 
I shall adopt here the bosonic picture. 

To emphasise the condensate, I write the grandcanonical partition function as a 
sum of terms corresponding to different values of no (the population of the condensate) 

n a n z—1 n D 

(9) 

The "partial" partition function, Z nQ , is proportional to the probability that a con- 
figuration with no particles in the condensate occures. As it happens usually at a 
first order phase transition, the partition function in some parameter space develops 
two maxima. In each of the phases one maximum is dominant and the equilibrium 
satate of the system corresponds to it. As we change from one phase to the other 
by varying some of the thermodynamic parameters of the system (e.g. the tempera- 
ture or the chemical potential), the difference between the two maxima decreases and 
they become equal at the phase transition. Moving further into the second phase, 
the maximum corresponding to it becomes dominant. For our system, the parameter 
with which we describe the phase transition is the population of the condensate, no. 

Therefore, instead of calculating the whole partition function, I calculate Z na 
and analyse its maxima. It is easier to work with In Z na , which is 

lnZ„ = + /3 M n + lnZ ox (n , (3). (10) 

a 

I introduce the dimensionless quantities a = (3/i, £ = f3N/a, and £o = Pno/a. In these 
notations Eq. (fTU]) becomes 

pa^taZn, =-h tl + ae, +pa- 1 \nZ ex {n ,p). (11) 

According to the calculation procedure outlined in Section [T] 2 ex (no,/3) is the par- 
tition function of a Fermi gas with the same chemical potential and with a single 
particle spectrum that starts at t\\ 

poo 

In Z cx (n , (3) = / In (l + e a ~^) ade. (12) 
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In Ref. [5] I showed that condensation may occure only on the first energy level, so I 
can write ?i = e% + wqct -1 = + no), which leads to 



/3a- 1 In Z cx (n , (3) = / In (l + e 

Jcr-^l+no) 

I plug this into (jTTj) and I get 



lnfl 



dx . 



(13) 



/•OO 

pa- 1 ]*^ = -h $+a£ + / ln^l + e"-^-^" 1 - 35 ) 



-ZioCo + a £o - Li 2 (-e 



), 



where I used the notation 



Li n (z) 



r(n) J e x z — 1 



(14) 



(15) 



for the polylogarithmic function of order n [151 116j . 

Now we have to find and compare the maxima of In Z ng . If In Z no forms two 
maxima, then our system has two phases. In each of the phases one maximum is 
dominant and no, together with all the measurable thermodynamic quantities, take 
the values corresponding to this maximum [17l HH [HI [20] ■ Since In Z no is a smooth 
function, we find its extrema at no > by equating to zero its derivatives with respect 
to no or, more conveniently, with respect to £o : 







dtio 



-2/ioCo + a - In 1 



We transform it a bit, to obtain 

/3a- 1 + toil - 2h ) =ln 



(16) 



(17) 



For an ideal Fermi gas, ho — 1/2 and we reobtain the equation for the Fermi conden- 
sate in ideal gases of constant a [T3] : 



Pa- 1 = In 



1 _|_ e -(a-£o)+/3c 1 



(18) 



In general, if ho > 1/2, the derivative d(Pa~ 1 In Z nQ ) / 'd^o is always negative and Z nQ 
has only one maximum at no = 0. We have therefore only one (uncondensed) phase 
for any a [S]- 

If ho < 1/2 we have a condensation 8 . In this case, since the condensate 
contains a macroscopic number of particles, we can neglect P/a as compared to a and 
£o and write Eq. (|17[) in the shorter form: 



£o(l-2fto) = ln[l 



(19) 



Equation (|19p is easy to interpret. For small enough a, the functions £o(l — 2/io) 
and ln[l + exp(£o — a)] do not cross each-other, so we have no solution for (|19p and 
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Figure 1: The partition functions of the canonical and grandcanonical ensembles and 
their derivatives 



Z no has only one maximum at uq = -no condensate. This situation is plotted with 
continuous lines in the right upper and lower plots of Fig. [TJ 

As a increases, Co(l — 2/irj) becomes tangent to ln[l + exp(Co ~ &)] when 

4- In[l + e 5 °- a l = = l-2h , 

which happends at Co,* = a — \n[2h n /(l — 2ho)]. Plugging Co,*, into Eq. (|T5|) I get 

2h ln(2/i ) + (1 - 2h ) ln(l - 2h ) 

at = T^2h • (20) 

This situation is drawn in dashed lines in the right-hand plots of Fig. [T] 

For a > a t , Eq. (|19|) has two solutions, Co,i and Co,2j with Co,i < £0,2- This 
means that for Co < £0,1 the partition function Z no decreases with increasing no, for 
£0,1 < Co < Co, 2, 2„ increases with no, and finally for Co > Co, 2, Z no decreases again. 
The two competing maxima of Z no , one at no = and the other one at no = fcBTcrCo.2, 
indicate a potential phase transition. Indeed, as a increases past at, the relative 
height of the maximum located at Co = Co, 2 increases with respect to the maximum 
at Co = 0, and at a certain critical value, a cr , we obtain Z no= o — Zn =k B Ta£ 1- 
Above a cr the second maximum is higher (see the dotted lines in the plots of Fig. [TJ , 
Z no= o = Z no= k B Tcr£o 2i an d we sa Y that a phase-transition occured. As we shall see 
in Section [5TT1 the (average) internal energy of the gas is conserved at the transition, 
but the particle number has a jump. 

This is a first order phase transition. 
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2.2 THE CANONICAL ENSEMBLE 

Let us now analyse the system in canonical conditions. Neglecting again the mi- 
croscopic energy step cr _1 in the system with interaction, the total particle number 
is 



N = n + N cx — n Q + 
If I denote £ G x = @<J~ l N cx , then 



/ aZa^-i ~n Q + k B Ta\n(l + e a -^) (21) 



C cx = ln(l + e Q -«°) (22) 

and Eq. (j2"Tj) becomes 

£ = £o + U = £o + In (1 + e Q -«°) . (23) 
In the canonical ensemble I express a in terms of £ and £oj 

a = £o + In (e*- & - 1) , (24) 
and then I write the parition function f3a~ 1 Z 7lQ in terms of £ and £o : 

poo 

$o- x Z na = (1 - ^oKo + Co In (e«-«« - l) + / In [l + (e«"«° - l) e"*] 

Jo 

= (1 - h )$ + £o In (e«-«« - l) - Li 2 (l - e«-«°) (25) 

Using Eq. I|25p I can find the maximum of j3a~ 1 Z no at fixed £, and for this I 
equate the derivative d([3a~ 1 Z no )/d£,o to zero, 

— ((3a- x Z no ) = (1 - 2/ioKo + In (l - e^"*) - £ ^ = . (26a) 

To discuss the solutions of Eq. I|26p , it is better to write it in the form 

(1 - 2fo )£o = _ 1 - In (1 - e*>-*) (26b) 

and denote the r.h.s. by /(£, £o)- The first and the second derivative of /(£, £o)> 

df(£, So) = + - 1 

3£o (e«-«o - l) 2 

and 



> 



^ ( e «-«o _ i) 3 



(^~ So + e s ~ £o + £ - 1) > 



are both bigger than zero for £ > and < £o < £; therefore f(£,£o), as a function 
of £o, is concave upwards. 
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To see if we have any solutions, I evaluate d((3a l Z no )/d£o at £o = and £o = £, 
and I get 



and 



n =Q 



In (l - e~«) 



e« - 1 



< 0, V£ > 0, (27a) 



0£o 



no—N 



= (1 - 2fto)€ + I™ 
= -oo, V£ > 0. 



(27b) 



In other words, the concave up-wards function /(£,£o) is bigger than the linear func- 
tion (1 — 2/io)£o a t both ends of the interval £o G [0>£]> so j if they are not tangent, 
they must cross each-other either zero times or two times in this interval. Therefore 
the Eqs. (|2"Sf have either no solutions or have two solutions for £ S (0, £). 

We have here the same situation as in the grandcanonical ensemble. For small 
£, Eq. (|26b[) have no solution. A situation like this is depicted in the left-hand 



plots in Fig. [T] with solid lines. As £ increases, at a certain value £t, the function 
(1 — 2/io)£o becomes tangent to /(£, £o), and -Z„ , as a function of £o, has an inflection 
at the tangent point, as shown with dashed lines in Fig. [T] For £ > £ t , Z„ has two 
maxima, one at £o = and another one at £o = £0,2 (I use the same notations as in 
the canonical ensemble). As £ increases above a critical value, say £ cr , the maximum 
formed at £0,2 becomes dominant (the dotted lines in Fig. [lj and at £ = £ cr we have 
a phase transition. 

In Fig. [2] I plot for comparison £0,2 as a function of /io for the grandcanonical 
(solid line) and canonical (dashed line) ensembles. It is remarcable that although 
as ho — ► 0.5 the Fermi gas approaches the ideal behavior, £0,2 does not converge to 
zero, as one would expect (i.e. no condensation in ideal, continous systems). This 
means that asymptotically, in the noninteracting system a finite fraction of the gas 
condenses at the transition temperature [l4] . 

In Section [3721 1 shall calculate the latent heat of this transition. 



3 DISCONTINUITIES AT THE PHASE TRAN- 
SITION 

Let us now calculate the jumps of the thermodynamic quantities that occure at the 
phase transition. Usually one expects that a first order phase transition in the grand- 
canonical enemble is marked by discontinuities in the internal energy (i.e. a latent 
heat) and in the particle number. In the canonical ensemble, the particle number is 
held fixes and the only quantity that may vary is the internal energy. If we think 
of typical first order phase transitions, like liquid- vapour transition or Bosc-Einstein 
condensation, we observe that they do not show any discontinuity in the internal 
energy in canonical conditions, but condense gradually as the temperature decreases. 
Unlike these, the transition described here, under canonical conditions, is marked by 
a jump of particles in the condensate, accompanied by another jump in the specific 
heat. 

I shall explain these in the next two subsections. 
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Figure 2: Comparison between of the jumps of the number of particles in the con- 
densate at the transition temperature, in canonical (dashed line) and grandcanonical 
(solid line) ensembles. 



3.1 GRANDCANONICAL ENSEMBLE 

Let us evaluate first the latent heat of the transition in the grandcanonical ensemble. 
If we denote by U the internal energy of the system and introduce the adimensional 
quantity U(a,£ l o) = f3 2 a~ 1 U, by standard calculations I obtain 

/>oo 

W(o, Co) = ho$ + £o In (1 + e°-«°) + / In (l + e -6-*) dx 

Jo 

= /i ^+^oln(l + e -f°) -Li 2 {-e a -t°) (28) 
Substracting l[14p form l[28p and ignoring /3<J _1 , I get 

Uia^-po-^hvZ^ = £ [ln(l + e«°- a ) - (1 - 2 V)£o] 

= Wo ' (29) 

so at the phase transition, U(a, £o) = Per" 1 In Z no , both at £o = and £o = Co, 2 (where 
Co, 2 corresponds again to the second maximum of ]xiZ no ). But since lniJ rao (a,£;o = 
0) = lnZ„ (a,Co,2) exactly at the phase transition, it implies that U(a, 0) = W(a,£o,2) 
and the grandcanonical latent heat, A 3 = (fcsT) 2 (r[iY(a, £0,2) — W(o, 0)] is identically 
zero. 

The quantity that varies at the grandcanonical phase transition is the particle 
number. From Eq. (j2"3")l I get directly 



S&ransition = £(a, ^0,2 ) ~ Co = 0) = In 



1 



1 + e- a 



> . (30) 
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Figure 3: The discontinuity of the particle number in the grandcanonical system at 
the transition temperature (continuous line). With dotted line is the 



The function S^transition is plotted in Fig. [3l together with the jump of the number 
of particles in the condensate, £0,2, as a function of ho. Unlike £0, £ converges to zero 
as ho approaches 0.5. 



3.2 CANONICAL ENSEMBLE 



To calculate the latent heat in the canonical ensemble, I replace a in Eq. 
the expression In this way I get 



/>oo 

U(Uo) = 6>(1 -/*>)]+ / In [1 + (e*-*° - 1) e" 

Jo 

= £o[£-£o(l- V)]-£*2(l-e«-«°) 
If I denote by A c the latent heat in the canonical ensemble, then 



dx 



-A c = £ [£-£o(i-M] 



In 



1 + (et-t° - l) 
1 + (e« - l)e 



dx . 



with 



(31) 



(32) 



This adimensional latent heat is plotted in Fig. [U 

Therefore the phase transition in the canonical system is marked by a disconti- 
nuity of the internal energy and a sudden condensation of a finite number of particles. 



4 CONCLUSIONS 



I analysed an interacting system model of constant density of single particle states, 
which can describe systems of generalised exclusion statistics. I showed that a first 
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0.1 0.2 0.3 0.4 0.5 

h 

Figure 4: The latent heat of the canonical ensemble, in units of (/cbT) 2 (7, as a function 
of the interaction parameter ho. 

order phase transition takes place and this phase transition is the condensation of a 
finite number of particles on the lowest energy levels. I characterised this transition 
for the situation when the system can be described as a Fermi gas. In grandcanonical 
conditions, the phase transition has no latent heat, but has a discontinuity in the 
number of particles, due to the condensation phenomenon. Unlike usual first order 
phase transitions, in this model we have a latent heat at the phase transition in 
canonical ensemble. 

Another interesting fact is that at the transition temperature, the number of 
particles in the condensate converges to a finite, nonzero value, as the system is 
changed gradually into a noninteracting system. This can be put into correspondence 
with the condensation that occures in ideal Fermi systems [14] , 
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